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■ The effects of dissipation on the thermodynamic properties of nonlinear quantum 
systems are approached by the path-integral method in order to construct approx- 
imate classical-like formulas for evaluating thermal averages of thermodynamic 
quantities. Explicit calculations are presented for one-particle and many-body 
systems. The effects of the dissipation mechanism on the phase diagram of two- 

■ , dimensional Josephson arrays is discussed. 
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CJ . 1 Introduction 

The usefulness of the improvedS effective potential approach, @ has been proven 
^ by several applications to condensed matter systems. However, open systems 

. were not immediately treated and previousLstudies were confined to obtain a 

classical-like expression for the free energy. I3c 
. In fact, the effective-potential method, also called pure-quantum self-con- 

sistent harmonic apmjDximatian (PQSCHA) after its generalization to phase- 
space Hamiltonians,|j'EI is ablelj to give the density matrix of a nonlinear systern 
(<~^ ' interacting with a dissipation bath through the Caldeira-Leggett (CL) model. □ 

For a better understanding of the method let us first consider one single degree 
' of freedom. In the general case the CL model starts from the Hamiltonian 



H=^ + V{q) + i V + m^ojl - F^) " L (1) 



rj^ • 2m ' ' 2 ^ 

Q , where {p,q) and {pi,qi) are the momenta and coordinates of the system and 

O ■ of an environment (or bath) of harmonic oscillators. When Fi{q) = q the 
dissipation is said to be linear and we will restrict ourselves to this case in 

, , ■ the following. The bath coordinates can be integrated out exactly from the 

r> I corresponding path integral and the CL Euclidean action is obtained in the 

■ form: 







S[q{u)] = / du ■-^e{u)+V{q[u))\ + 1 - jdu' k{u~u') q[u) q{u') . (2) 
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To make contact with the classical concept of dissipation, one may take the 
classical counterpart of Eq. and get the Langevin equation of motion: 

/>oo 

mq + m dt' -f{t') q{t-t') + V'(q) ; 



then the Laplace transform 7(2) of the damping function j{t)-as expressed in 
terms of the spectral density of the environmental coupling, can be related 
to the Matsubara components of the damping kernel k{u) {vn = i'Kn/ (3) : 

k{u) = E„ e*'"" , = \l^n\l{z=Wn\) , (3) 

We will consider two cases: 

Ohmic : '-f{t) = j S{t — 0) , 7(2) — 7 , 

Drude : -/{t) ^ uju e^"^"*, 7(2) = 7 cjD/(t^D + 2) , 

where the dissipation strength 7 and the bath bandwidth lod characterize the 
environmental coupling; wd — * 00 gives the Ohmic (or Markovian) case. In 
order to interpolate between different regimes and to take into account memory 
effects we can also assume: 

7(z) oc z' , (4) 
with — 1 < s < 1 and s > (s < 0) is said the super- (sub-) Ohmic case. 

2 Effective potential in presence of dissipation 

The PQSCHA approximation consists in taking as trial action So, the most 
general quadratic functional with the same linear dissipation of S, namely, 

So[q{u)] ^ / du —q^{u) + w + m,Lo'^{q{u)~qy 
.In L ^ 



du 
T 



du' k{u—u') q{u) q(u') . (5) 



The quantity q = (3 ^ J q{u) du is the average point of paths and 

w = w{q) , UJ^ = uj^{q) , (6) 

are parameters to be determined by minimizing the right-hand side of the 
Feynman inequality, 

F<Fo+r'{S^So)so ■ (7) 
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For any obseiyable 0{p, q) the S'o-average can be expressed Q'Ll in terms of its 
Weyl symbol □ q). As final result we obtain the classical- like form 

<'*>^i\/S/'""«^(''-'+«»>^""""""- 

where ((•••)) is a Gaussian average operating over p and ^ with moments 

^ °fa-) = ^ g .g+„.(,-)H.t„ — 2^^(»"'^-7) ' O 
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where / = /3w/2. The effective potential is defined as Vcs{q) = w{q) + a{q) , 
with 



The r.h.s. of Eqs. (p|-|l 1|) show that the well known non-dissipative limits 
are recovered for fc — > 0. The variational parameters can be self-consistently 
calculated and the explicit expressions are the following 

w{q) = {{V{q + 0)) - \rnoj^a , muj^q) = {{V"{q + i))) . (12) 

3 Applications 

In order to understand how this approximation scheme works, let us take 
first a very simple system: one particle in a double-well potential with Ohmic 
dissipation. A typical result for this system is shown in Fig. |l]: the dissipation 
quenches the quantum fluctuations of the coordinate. However, it must be 
noted from Eq.(|l0|) that those of the momentum are infinite. 

Let us now turn to the many-body case. The fixs,t application is a quantum 
</?^-chain of particles with Drude-like dissipation: Ej the undamped system is 
described by the following action 



2a2 8 



(13) 



where a is the chain spacing, is the gap of the bare dispersion relation, 
and Q is the quantum coupling. The classical continuum model supports 
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kink excitations of characteristic width fl^^ and static energy ek — fi/Q. 
This energy is used as the energy scale in defining the reduced temperature 
t = {f3eK)~^- We also use the kink length in lattice units R = {a^l)~^ (i? — > c» 
in the continuum limit). We assume uncorrelated identical GL baths for each 
degree of freedom, and use the low-coupling approximation □ (LCA) in order 
to deal with the effective potential. The partition function turns out to be 
written as 



3t 



q e 
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(14) 
(15) 
(16) 



The renormalization parameter D = D{t;Q,R;^,ujY}) generalizes Eq. (^ 
and is the solution of the self-consistent equations: 



D 



ARt 
3iV 



E E - 



Lolit) = f^^ [(1 - iD) + 4i?2 sin2(fca/2)] , 
where fc„ — 70^0 i^nj ("-^D + i^n) and here z/„/r2 = lixntlQ. 



(17) 
(18) 




Figure 1: Configuration density P(x) = — x)) of the double-well quartic potential for 
the fixed eoupling g = 5, the redueed temperature < = 1, and different values of the Ohmic 
damping parameter F = 7/(^0, being ci^o the characteristic frequency of the system. The 
filled circles are the exact result for F = 0; the dotted curve at F = 00 corresponds to the 
classical limit. 
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Again, we observe that the fluctuations of coordinate-dependent observ- 
ables are quenched by increasing the dissipation strength 7, while those of 
momentum-dependent ones are enhanced due to the momentum exchanges 
with the environment. The result of these opposite behaviors is a non-trivial 
dependence on dissipation of "mixed" quantities like the specific heat as is 
shown in Fig. ^ . 

The last system we consider is the dissipative quantum XY model. Much 
interest in such system is due to its close relation with 2D granular super- 
conductors and Josephson junction arrays (JJA). In the last 15 years nuich 
attention has been devoted to the theoretical and experimental study O of 
the phase ordering in JJA and of how it is influenced by the environmental 
coupling. The undamped system is described by the action 

5jjA = Td^l^E'^^HTl'^^^") -^T. cos[Mu)-<Pj{u)]\ , (19) 
"'0 ij ^ <ij> ^ 

where (pi is the superconducting phase of the i-th island, <zj> restricts the 
sum over nearest-neighbor bonds, J is the Josephson coupling, dj = (Co -I- 

ACi)5ij — Ci(5f^"\ Co and Ci are the self and mutual capacitance between 

superconducting islands and 5^™^ — 1 for nearest neighbors, zero otherwise. 

From Eq. (|l9|) two contributions to the energy in this system can be ob- 
served: the kinetic one, due to the charging energy between islands [E^ = 
q = 2e), and the potential energy, due to the Josephson coupling in the 




Figure 2: Total specific heat c{t) and the kinetic and interaction parts of the specific heat, 
namely c^, (t) = dt (K) and (t) = dt {V) vs. reduced temperature t, for different values of 
the damping strength F = 7/n. Note that for F — > 00, c{t) tends to c^i — 1/2: this behavior 
can be explained by observing that, in the strong damping limit, c^{t) 0. 
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superconducting junctions. When Ec<S^J, the charges on each island fluctuate 
independently from the phases (pi : the latter have then a classical XY behav- 
ior and the associated Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
takes place at temperature Tj^'^)^ ^ 0.892 J. In the opposite limit, the energy 
cost to transfer charges between neighboring islands is too high, so the charges 
tend to be localized and the phase ordering tends to be suppressed. 

In this scenario dissipation is supposed to have an important role^ The 
environmental interaction tends to suppress the quantum fluctuations E3 of (f>i 
and to restore an almost classical BKT phase transition. Nevertheless, it is 
not clear which is the physical mechanism of the dissipation. For the single 
Josephson junction with Ohmic dissipation the classical "resistively and capac- 
itively shunted junction" (RCSJ) model is recoveredtS, but in the case of many 
degrees of freedom the environmental interaction is much more complicated, 
e.g. non-exponential memory effects and then non-Ohmic damping can ap- 
pear. The dissipation model we assume consists in independent environmental 
baths, one for each junction (or bond). The dissipative part of the action is 

'^nKn4>n, (20) 

where the Fourier transform of the CL kernel matrix is given by 

and LOa is a characteristic frequency, that we choose as the Debye frequency 

/ ^2 J \ 1/2 (J 

Now it is possible to define the quantum coupling parameter, that measures the 
"quanticity" of the system as the ratio between the characteristic quantum and 
classical energy scales g — hujo/J. The effective potential is calculated using 
the extension to the many-body casellj of the scheme of Sec. ^and, apart from 
uniform terms, is given by 

Ves = - Jaff(r) J2 cos(0. -(pj) , (23) 
<y> 

where Joff(T) = Je"^^*^'^-'/^ and Di{T) is the renormalization parameter 
that measures the pure-quantum contribution to the fluctuations of the nn 
relative superconducting phases. The phase diagram of the system is cal- 
culated starting from the classical effective potential (^3|), i.e. by solving 
"^BKT — Jcs{T^ii^)T^^^ , and is shown in Fig. ||, for different values of the 
parameters that characterize the dissipation. 
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Figure 3: Phase diagram in the (<?,<) plane, t = T/J, for different values of the damping 
parameters, T = 7/0)0 and s. On the left, the case of Ohmic dissipation: the critical 
temperature T^j^.j, tends to T^kt ^^"^ increasing damping strength. On the right, for different 
values of s: the cases s = 1 and s = —I are nondissipative, since they correspond to a 
variation of the capacitance- and of the frequency spectrum, respectively. 



7 



